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PEPSIG  (CORNER  VORTEX  VERSION) 
INPUT  DESCRIPTION 


Input  to  PEPSIG  (Corner  Vortex  version)  consists  of  an  initial  card  specifying  the 
type  of  calculation  to  perform  and  a  title,  plus  three  namelists.  Many  of  the  input 
parameters  have  default  values  and  do  not  need  to  be  specified  by  the  user  unless 
some  other  value  is  desired.  The  type  (REAL  or  INTEGER)  of  the  input  parameters 
follows  standard  FORTRAN  convention,  unless  stated  otherwise.  (I.e.,  those  starting 
with  I,  J,  K,  L,  M,  or  N  are  INTEGER,  and  the  remainder  are  REAL). 


TITLE  CARD 

MODE  An  integer  in  column  1  specifying  the  type  of  calculation  to  be 

performed,  as  follows: 

3  to  compute  and  store  geometric  parameters  for  potential 

flow  calculation. 

2  to  perform  potential  flow  calculation. 

1  to  perform  viscous  flow  calculation. 

0  to  terminate  xun. 

TITLE  A  descriptive  title  in  columns  2-33  used  on  the  printed  output  and  in 

the  plot  file. 


NAMELIST  RESTRT 

The  parameters  specified  in  this  namelist  are  primarily  used  to  control  where  a 
restart  file  is  read  and/or  written.  A  restart  file  must  be  used  when  changing  the 
polynomials  used  to  specify  the  geometry,  the  marching  step  size  DT,  the  type  of 
output  desired,  etc. 

IRSTIN  Marching  station  number  to  be  read  from  restart  file.  IRSTIN  =  0 

implies  this  is  not  a  restart  case.  The  default  value  is  0. 

IRSTOT  The  interval  for  writing  onto  the  restart  file  (i.e.,  at  station 

IRSTIN+IRSTOT,  IRSTIN+2*IRSTOT,  etc.).  If  IRSTOT  =  0  no  restart  file 
is  written.  The  default  value  is  0. 

NFILE  The  sequence  number  in  the  restart  file  of  the  station  to  be  read. 

E.g.,  if  stations  1,  5,  20,  and  24  have  been  written  onto  the  restart 
file,  and  a  restart  at  station  20  is  desired,  then  NFILE  should  be  3. 
The  default  value  is  1. 

NSAVED  The  number  of  stations  in  the  restart  file  to  be  saved.  E.g.,  if 

stations  1,  5,  20,  and  24  have  been  written  onto  the  restart  file,  and 
NSAVED  =  3,  then  when  a  restart  is  next  written  stations  1,  5,  and  20 
will  be  saved  and  station  24  will  be  overwritten.  The  default  value  is 
NFILE.  Note  that  if  JRSTIN  =  JRSTOT  (see  below)  and  NFILE  =  NSAVED, 
the  same  file  can  be  used  for  reading  and  writing  restart  data  with 
destroying  any  previously  saved  data. 


NAVY  USERS  MANUAL 


INPUT 


2 


NOTE:  In  practice  NFILE  and  NSAVED  are  both  usually  defaulted, 
which  results  in  each  write  to  the  restart  file  overwriting  the 
previously  saved  station. 

JRSTIN  Fortran  unit  number  from  which  restart  data  will  be  read.  The 

default  value  is  11. 

JRSTOT  Fortran  unit  number  onto  which  restart  data  will  be  written.  The 

default  value  is  11. 


NAMELIST  FLUIDS 

The  first  group  of  parameters  in  this  namelist  apply  to  all  cases,  and  are  used  to  set 
reference  conditions,  initial  profiles,  and  boundary  conditions. 

IUNITS  1  Dimensionless  input  and  output. 

2  SI  units  in  input  and  output. 

3  English  units  in  input  and  output. 

The  default  value  is  1. 

INPOPT  Eight  options  are  available  for  specifying  a  consistent  set  of 
reference  conditions.  These  conditions  are  also  used  to  set  up  the 
initial  profiles  to  start  the  viscous  marching  calculation.  For 
INPOPT  =  1  to  8,  the  parameters  to  be  specified  by  the  user  are 
summarized  as  follows: 

NPOPT  PARAMETERS  SPECIFIED 

1  CMACH  REY  YZERO  TZERO  PZERO 

2  UZERO  REY  YZERO  TZERO  PZERO 

3  CMACH  REY  YZERO  TZERO  RZERO 

4  UZERO  REY  YZERO  TZERO  RZERO 

5  CMACH  REY  YZERO  SOUND  PZERO 

6  UZERO  REY  YZERO  SOUND  PZERO 

7  CMACH  REY  YZERO  SOUND  RZERO 

8  UZERO  REY  YZERO  SOUND  RZERO 

The  default  value  is  1. 

KTURB  0  Laminar  flow. 

1  Turbulent  flow. 

The  default  value  is  0. 

CMACH  Reference  Mach  number.  The  default  value  is  0.01. 

REY  Reference  Reynolds  number,  based  on  RZERO,  UZERO,  and 

YZERO.  The  default  value  is  1. 

UZERO  Reference  velocity.  The  default  value  is  1. 

YZERO  Reference  length.  The  default  value  is  1. 


□  □ 
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TZERO  Reference  static  temperature.  The  default  value  is  1. 

PZERO  Reference  static  pressure.  The  default  value  is  1. 

RZERO  Reference  static  density.  The  default  value  is  1. 

SOUND  Reference  speed  of  sound.  There  is  no  default  value. 

RG  Gas  constant.  The  default  value  is  1716.,  which  works  for  both 

IUNITS  =  1  and  3. 

CPR  Reference  value  for  specific  heat  at  constant  pressure,  used  in 

formula  for  specific  heat  as  a  function  of  temperature.  The  default 
value  is  6006.,  which  also  works  for  both  IUNITS  =  1  and  3. 

BLD  A  6-element  array  giving  the  initial  boundary  layer  thicknesses  on 

the  four  computational  surfaces,  and  on  the  inner  and  outer  lobe 
surfaces  for  a  mixer  case.  The  default  values  are  0.,  0.,  0.,  0.2,  0.,  and 
0. 

VGEN(J,20)  Coefficients  of  the  polynomial  specifying  the  boundary  layer 
thickness  along  the  vertical  wall.  The  independent  variable  in  this 
polynomial  is  the  height,  z.  Up  to  a  second-order  polynomial  may 
be  specified  (i.e.,  J  =  1  to  3).  J  =  n  corresponds  to  the  z  to  the  (n-l)’th 
term  in  the  polynomial. 

VGEN(J,19)  Coefficients  of  the  polynomial  specifying  the  boundary  layer 
thickness  along  the  horizontal  wall.  The  independent  variable  in 
this  polynomial  is  the  distance  from  the  corner,  y.  Up  to  a  second- 
order  polynomial  may  be  specified  (i.e.,  J  =  I  to  3).  J  =  n  corresponds 
to  the  y  to  the  (n-l)'th  term  in  the  polynomial. 

QRATIO  For  internal  flows  -  ratio  of  the  average  initial  velocity  at  the  start 
of  the  viscous  run  to  the  inlet  velocity  in  the  potential  flow  run. 
This  parameter  is  used  in  setting  the  potential  flow  pressure  field, 
and  is  most  important  when  the  initial  station  in  the  viscous  run 
does  not  have  the  same  mass  flow  as  the  initial  station  in  the 
potential  flow  run.  To  determine  what  value  to  use:  (1)  set  QRATIO  = 
1.0  and  run  the  viscous  calculation  one  step;  (2)  check  the  printed 
output  for  the  value  of  the  average  velocity  at  the  initial  station 
(UBAR);  (3)  run  the  complete  viscous  case  with  QRAT  -  UBAR.  The 
default  value  is  1. 


NS1 

0 

Symmetry 

boundary 

condition 

on 

computational 

surface 

1 

1 

Solid  wall 

boundary 

condition 

on 

computational 

surface 

1. 

The 

default  value  is  0. 

NS2 

0 

Symmetry 

boundary 

condition 

on 

computational 

surface 

2 

1 

Solid  wall 

boundary 

condition 

on 

computational 

surface 

2. 

The 

default  value  is  0. 

NS3 

0 

Symmetry 

boundary 

condition 

on 

computational 

surface 

3 

1 

Solid  wall 

boundary 

condition 

on 

computational 

surface 

3. 
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The  default  value  is  0. 

NS4  0  Symmetry  boundary  condition  on  computational  surface  4. 

1  Solid  wall  boundary  condition  on  computational  surface  4. 

The  default  value  is  1. 

ISYM  Degree  of  symmetry  in  the  cross-section,  i.e.,  ISYM  =  1  for  a  full  360 

degree  polar  cross  section,  ISYM  =  2  for  a  90  degree 

calculation, ISYM  =  4  for  a  90  degree  calculation,  etc.  The  default 
value  is  2. 

IPRD  0  No  periodic  boundary  condition. 

1  Periodic  boundary  condition  in  the  circumferential  direction. 
The  default  value  is  0. 

NOTE  -  IPRD  =  1  only  applies  to  the  NGEOM  =  2-5  and  12-15  geometry 
options.  In  addition,  when  IPRD  =  1  stretched  meshes  in  the 
circumferential  direction  cannot  be  used. 

ILAW  0  Compute  the  laminar  viscosity  coefficient  from  Sutherland's 

law. 

1  Compute  the  laminar  viscosity  coefficient  by  assuming  it  is 
proportional  to  temperature  to  the  0.76  power. 

2  Hold  the  laminar  viscosity  coefficient  constant  at  its  input 
value  (computed  from  YZERO,  UZERO,  RZERO,  and  REY). 

The  default  value  is  2. 

IVT  0  Hold  the  specific  heat  at  constant  pressure  (CP)  and  the  ratio  of 

specific  heats  (GAMMA)  constant  at  their  input  values 
(computed  from  TZERO). 

1  Treat  CP  and  GAMMA  as  variable  functions  of  temperature. 

The  default  value  is  0. 

FLRFC  Flare  approximation  criteria  for  separated  flow  regions.  If  U  < 

FLRFC*UAVG,  U  is  reset  to  FLRFC*UAVG.  The  default  value  is  0.025. 

The  following  parameters  also  apply  to  all  cases,  and  are  used  to  specify  the 

distribution  of  grid  points  in  the  streamwise  direction,  and  the  number  of  grid 

points  in  the  transverse  directions.  The  distribution  of  grid  points  in  the 

transverse  directions  is  controlled  by  parameters  in  namelist  GEOM. 

T  Marching  parameter,  or  streamwise  computational  coordinate. 

This  is  the  independent  variable  used  in  evaluating  the 
polynomials  PCL  and  PGEO  in  namelist  GEOM.  The  value  of  T  must  be 
specified  at  the  initial  station.  After  a  restart,  however,  the  value  is 
taken  from  the  restart  file.  There  is  no  default  value. 

DTE(K)  Step  size  in  T  for  marching  step  from  station  1RSTIN  to  IRSTIN+1. 

There  is  no  default  value.  Up  to  10  piecewise  continuous  sections, 

each  beginning  at  TSECT(K),  can  be  specified,  (default  K=l). 

AP(K)  Ratio  of  successive  step  sizes  in  T.  I.e.,  AP  =(T(I  +  1)-T(I))/(T(I)-T(I- 

1)).  The  default  value  is  1.0,  corresponding  to  a  constant  step  size. 
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Up  to  10  piecewise  continuous  sections,  each  beginning  at  TSECT(K), 
can  be  specified,  (default  K=l). 

NS  Number  of  last  streamwise  station  to  be  computed.  I.e.,  the  code  will 

march  from  station  IRST1N  to  station  NS.  There  is  no  default  value. 

NEY  Number  of  grid  points  in  the  circumferential  direction.  There  is  no 

default  value. 

NEZ  Number  of  grid  points  in  the  radial  direction.  There  is  no  default 

value. 

The  following  parameters  apply  to  all  cases  and  control  the  type  and  amount  of 
output  to  be  printed  and/or  written  to  the  plot  file.  (ICOEF  also  controls  other 
aspects  of  the  calculation). 

IOOEF  A  6x20  array  containg  various  switches,  limits,  etc.  See  the  separate 

write-up  on  ICOEF  for  the  details. 

IWSTA  A  101-element  array  specifying  station  numbers  at  which  output 

will  be  printed.  This  option  is  activated  when  1C0EF(2,3)  <  0.  The 
default  values  are  all  0. 

KPRT(I)  An  array  of  up  to  15  elements  used  to  specify  additional  printout 

beyond  the  standard  printout.  This  is  done  by  setting  KPRT(I)  equal 
to  an  integer  corresponding  to  the  variable  desired.  For  values 

from  1  to  MV  ARP  (a  PARAMETER  equal  to  the  number  of  variables 
in  the  Z  array),  the  corresponding  variable  from  the  Z  array  is 
printed.  Additional  variables  may  also  be  printed.  Table  ?  lists  the 
variables  stored  in  the  Z  array,  plus  the  others  that  may  be  printed. 
For  example,  setting  KPRT(l)  =  62  and  KPRT(2)  =  71  results  in 
profiles  of  turbulence  mixing  length  and  Mach  number  being 

printed. 

The  default  values  are  all  0. 

IPLOT  0  Do  not  write  a  plot  file. 

1  Write  a  file  for  later  post-processing.  The  stations  written  into 
the  file  are  controlled  by  the  values  of  ICOEF(2,4),  ICOEF(2,5), 
and  IPSTA. 

The  default  value  is  1. 

IPSTA  A  101 -element  array  specifying  station  numbers  at  which  output 

will  be  written  into  the  plot  file.  This  option  is  activated  when 

ICOEF(2,5)  <  0.  The  default  values  are  all  0. 

NPLT  Number  of  variables  written  into  the  plot  file  at  each  station.  The 

maximum  allowed  is  MVARP  +  10.  The  default  value  is  14. 

KPLT(I)  An  array  of  NPLT  elements  specifying  which  variables  are  to  be 
written  into  the  plot  file.  This  is  done  in  the  same  way  as  the 
specification  of  additional  printout  using  KPRT.  The  defaults  for 
KPLT  are  such  that  the  following  14  variables  are  stored  in  the  plot 
file: 
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1  Reference  Cartesian  coordinate  in  x-direction.  Here 

"reference  coordinate"  refers  to  a  fixed  coordinate  system,  as 
opposed  to  one  that  is  perpendicular  to  a  duct  centerline,  for 
example. 

2  Reference  Cartesian  coordinate  in  y-direction. 

3  Reference  Cartesian  coordinate  in  z-direction. 

4  Inviscid  velocity  (from  potential  flow  file)  in  relative 

Cartesian  y-direction. 

5  Inviscid  velocity  (from  potential  flow  file)  in  relative 

Cartesian  z-direction. 

6  Static  pressure  from  solution  of  Poisson  pressure  equation. 

7  Static  pressure  coefficient. 

8  Inviscid  static  pressure  coefficient  (from  potential  flow  file). 

9  Streamwise  gradient  of  inviscid  static  pressure. 

10  Velocity  in  relative  Cartesian  y-direction,  non-dimensionalized 

by  UZERO.  Here  "relative  Cartesian"  refers  to  a  coordinate 
system  whose  x-axis  is  normal  to  the  plane  of  the 
computational  station  (and  therefore  the  y-z  axes  are  in  the 
plane).  The  x-direction  is  often  referred  to  as  the  streamwise 

or  primary  flow  direction. 

1 1  Velocity  in  relative  Cartesian  z-direction,  non-dimensionalized 

by  UZERO. 

12  Velocity  in  relative  Cartesian  x-direction  (i.e.,  streamwise 

direction),  non-dimensionalized  by  UZERO. 

13  Vorticity  component  in  the  streamwise  direction,  non- 

dimensionalized  by  UZERO/YZERO. 

14  Secondary  stream  function,  non-dimensionalized  by 

UZERO*YZERO. 

The  following  parameters  apply  to  all  cases,  and  control  the  starting  procedure  for 
the  initial  marching  step,  and  the  sequence  of  solution  of  the  equations. 

KSTART  0  Original  PEPSIG  starting  procedure. 

1  Iterative  starting  procedure. 

2  Split  the  first  marching  step  into  substeps. 

The  default  value  is  0. 

DXSTRT  Initial  marching  parameter  step  size  for  KSTART  =  1  option.  The 

default  value  is  T+AP*DT. 
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NSTART  Number  of  iterations  in  the  starting  sequence  if  KSTART  =  1,  or 
number  of  substeps  if  KSTART  =  2. 

The  default  value  is  1. 

1SEQ  0  Original  equation  solving  sequence  in  PEPSIG. 

1  Optional  equation  solving  sequence. 

2  Another  optional  equation  solving  sequence. 

The  default  value  is  0. 

IBETA  0  Use  the  small  scalar  approximation. 

1  Use  the  pressure  approximation. 

The  default  value  is  1. 


The  remaining  variables  in  namelist  FLUIDS  are  either  not  used,  apply  to  options 
currently  being  developed,  or  are  for  output  only. 


NAMELIST  GEOM 


The  following  parameters  specify  the  geometry  being  analyzed. 


NGEOM 


2  Superelliptic  cross-sections,  set  up  through  subroutine  TUBE, 
with  the  semi-major  axis,  ratio  of  semi-major  to  semi-minor 
axis,  and  exponent  specified. 

3  Superelliptic  cross-sections,  set  up  through  subroutine  TUBE, 
with  the  semi-major  axis,  semi-minor  axis,  and  exponent 
specified. 

4  Same  as  NGEOM  =  2  but  with  elliptical  centerbody. 

5  Same  as  NGEOM  =  3  but  with  elliptical  centerbody. 

12  Same  as  NGEOM  =  2  but  going  through  ANYGEO  and  GEOTUB 
instead  of  TUBE. 

13  Same  as  NGEOM  =  3  but  going  through  ANYGEO  and  GEOTUB 
instead  of  TUBE. 

14  Same  as  NGEOM  =  4  but  going  through  ANYGEO  and  GEOTUB 
instead  of  TUBE. 

15  Same  as  NGEOM  =  5  but  going  through  ANYGEO  and  GEOTUB 
instead  of  TUBE. 

21  Rectangular  cross-sections. 

The  default  value  is  2. 


NOTE:  The  NGEOM  =  2-5  options  will  eventually  be  replaced  by  the  NGEOM  =  12-15 
options.  Similarly,  the  NGEOM  =  21  option  may  also  eventually  be 
replaced  by  an  ANYGEO  equivalent.  See  the  separate  write-up  on 
ANYGEO  for  details  on  how  to  add  additional  user-written  geometry 
packages. 

TSECT(K)  The  geometry  can  be  specified  in  up  to  10  piecewise  continuous 

sections.  Section  K  begins  at  TSECT(K)  (default  K=l). 

PCLD(J,1,K)  Coefficients  of  the  polynomial  specifying  the  y  Cartesian 

component  of  the  centerline  location,  if  ICOEF(l,9)  =  1;  the 

curvature  CNU  of  the  centerline  (i.e.,  the  reciprocal  of  the  local 
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radius  of  curvature),  if  IC0EF(1,9)  =  0;  or  the  Frenet  frame 
curvature  CAPPA  of  the  centerline,  if  ICOEF(l,9)  ■=  2.  The 
independent  variable  in  this  polynomial  is  the  marching 
parameter  T.  Up  to  a  tenth-order  polynomial  may  be  specified 
(i.e.,  J  =  1  to  11).  J  =  n  corresponds  to  the  x  to  the  (n-l)'th  term  in 
the  polynomial.  The  default  values  are  all  0,  corresponding  to  a 
straight  centerline.  Up  to  10  piecewise  continuous  sections,  each 

beginning  at  TSECT(K),  can  be  specified,  (default  K=l). 

PCLD(J,2,K)  Coefficients  of  the  polynomial  specifying  the  z  Cartesian 

component  of  the  centerline  location,  if  1C0EF(1,9)  =  1;  the 

curvature  CETA  of  the  centerline  (i.e.,  the  reciprocal  of  the  local 
radius  of  curvature),  if  ICOEF(l,9)  =  0;  or  the  Frenet  frame  torsion 
TAU  of  the  centerline,  if  ICOEF(l,9)  =  2.  The  independent  variable 
in  this  polynomial  is  the  marching  parameter  T.  Up  to  a  tenth- 
order  polynomial  may  be  specified  (i.e.,  J  =  1  to  11).  J  =  n 

corresponds  to  the  x  to  the  (n-l)’th  term  in  the  polynomial.  The 
default  values  are  all  0,  corresponding  to  a  straight  centerline.  Up 
to  10  piecewise  continuous  sections,  each  beginning  at  TSECT(K), 

can  be  specified,  (default  K=l). 

PCLD(J,3,K)  Coefficients  of  the  polynomial  specifying  the  x  Cartesian 

component  of  the  centerline  location,  if  ICOEF(l,9)  =  1,  or  the  arc 

length  along  the  centerline,  if  ICOEF(l,9)  =  0  or  2.  The 
independent  variable  in  this  polynomial  is  the  marching 
parameter  T.  Up  to  a  tenth-order  polynomial  may  be  specified 
(i.e.,  J  =  1  to  11).  J  =  n  corresponds  to  the  x  to  the  (n-l)’th  term  in 
the  polynomial.  The  default  values  are  all  0,  except  PCL(2,3)  =  1., 
which  says  that  the  marching  parameter  T  is  the  Cartesian  x 
coordinate.  Up  to  10  piecewise  continuous  sections,  each 
beginning  at  TSECT(K),  can  be  specified,  (default  K=l). 

PGEOD(J,l,K)  Coefficients  of  the  polynomial  specifying  the  ratio  of  the  semi¬ 

major  to  semi-minor  axis,  if  NGEOM  =  2,  4,  32,  or  34,  or  the  semi¬ 
minor  axis  itself  if  NGEOM  =  3,  5,  33,  or  35.  The  independent 
variable  in  this  polynomial  is  the  marching  parameter  T.  Up  to  a 
tenth-order  polynomial  may  be  specified  (i.e.,  J  =  1  to  11).  J  =  n 
corresponds  to  the  x  to  the  (n-l)'th  term  in  the  polynomial.  The 
default  values  are  all  0,  except  PGEO(l,l)  =  1.  Up  to  10  piecewise 
continuous  sections,  each  beginning  at  TSECT(K),  can  be  specified, 
(default  K=l). 

PGEOD(J,2,K)  Coefficients  of  the  polynomial  specifying  the  semi-major  axis. 

The  independent  variable  in  this  polynomial  is  the  marching 
parameter  T.  Up  to  a  tenth-order  polynomial  may  be  specified 
(i.e.,  J  =  1  to  11).  J  =  n  corresponds  to  the  x  to  the  (n-l)’th  term  in 
the  polynomial.  The  default  values  are  all  0,  except  PGEO(l,2)  =  1. 
Up  to  10  piecewise  continuous  sections,  each  beginning  at 
TSECT(K),  can  be  specified,  (default  K=l). 

PGEOD(J,4,K)  Coefficients  of  the  polynomial  specifying  the  exponent  of  the 

superellipse.  The  exponent  must  be  greater  than  or  equal  to  2.0 
and  should  be  less  than  or  equal  to  10.0.  The  independent  variable 
in  this  polynomial  is  the  marching  parameter  T.  Up  to  a  tenth- 
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order  polynomial  may  be  specified  (i.e.,  J  =  1  to  11).  J  -  n 
corresponds  to  the  x  to  the  (n-l)'th  term  in  the  polynomial.  The 

default  values  are  all  0,  except  PGEO(l,4)  =  1.  Up  to  10  piecewise 
continuous  sections,  each  beginning  at  TSECT(K),  can  be  specified, 

(default  K=l). 

PGEOD(J,5,K)  Coefficients  of  the  polynomial  specifying  the  semi-major  axis  of 

the  elliptical  centerbody.  The  independent  variable  in  this 
polynomial  is  the  inarching  parameter  T.  Up  to  a  tenth-order 

polynomial  may  be  specified  (i.e.,  J  =  1  to  11).  J  =  n  corresponds  to 
the  x  to  the  (n-l)'th  term  in  the  polynomial.  The  default  values 
are  all  0.  Up  to  10  piecewise  continuous  sections,  each  beginning 

at  TSECT(K),  can  be  specified,  (default  K=l). 

PGEOD(J,6,K)  Coefficients  of  the  polynomial  specifying  the  semi-minor  axis  of 

the  elliptical  centerbody.  The  independent  variable  in  this 
polynomial  is  the  marching  parameter  T.  Up  to  a  tenth-order 

polynomial  may  be  specified  (i.e.,  J  =  1  to  11).  J  =  n  corresponds  to 
the  x  to  the  (n-l)’th  term  in  the  polynomial.  The  default  values 
are  all  0.  Up  to  10  piecewise  continuous  sections,  each  beginning 

at  TSECT(K),  can  be  specified,  (default  K=l). 

Radius  of  small  tube  surrounding  polar  coordinate  singularity  at 
the  centerline.  The  recommended  value  is  half  the  distance  to  the 
first  grid  point  from  the  centerline.  The  default  value  is  0.01. 

Centerline  arc  length  at  the  start  of  the  calculation. 

The  default  value  is  0.0. 

0  External  flow  case.  (Requires  input  reference  point  for 
pressure.  See  ICOEF(2,16)  and  (2,17)). 

1  Internal  flow  case. 

The  default  value  is  1. 

0  Print  grid  point  locations  in  relative  Cartesian  coordinates. 

1  Print  grid  point  locations  in  polar  coordinates. 

The  default  value  is  1. 

ITUBE  is  automatically  set  equal  to  0  if  NGEOM  >  20. 

The  following  parameters  control  the  distribution  of  grid  points  in  the  transverse 
directions. 

IPA 


1PB 


1  Pack  grid  points  in  the  circumferential  direction  using 

subroutine  KGRID,  as  specified  by  NCTRK,  XAK,  and  EPSK. 

2  Uniformly  spaced  circumferential  grid. 

3  Pack  grid  points  in  the  circumferential  direction 

automatically,  based  on  the  cross-section  geometry. 

The  default  value  is  2. 

1  Pack  grid  points  in  the  radial  direction  using  subroutine 

KGRID,  as  specified  by  NCTRK,  XAK,  and  EPSK. 

2  Uniformly  spaced  radial  grid. 

3  Pack  grid  points  in  the  radial  direction  as  specified  by  VIS. 
The  default  value  is  2. 


RTUBE 

ARC 

IDUCT 

ITUBE 
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VIS  Radial  grid  packing  parameter  used  when  IPB  =  3.  A  value  of  0.15 

gives  a  loosely  packed  grid  near  the  outer  wall,  and  a  value  of  0.6 

gives  a  very  tightly  packed  grid  near  the  outer  wall.  For  laminar 

flow,  values  of  0.15  to  0.4  are  recommended.  For  turbulent  flow, 
values  of  0.4  to  0.6  are  recommended.  The  default  value  is  1.  (why, 
I  don't  know). 

NCTRK(J)  Number  of  locations  between  end  points,  in  transverse 

computational  coordinate  directions  (J  =  1  and  J  =  2)  about  which 

grid  points  are  to  be  packed  using  KGRID.  The  maximum  value 
allowed  is  10.  There  are  no  default  values. 

XAK(I,J)  Non-dimensional  locations  (i.e.,  from  0.0  to  1.0),  in  transverse 

computational  coordinate  directions  (J  =  1  and  J  =  2)  about  which 

grid  points  are  to  be  packed  using  KGRID.  XAK(1,J)  and 
XAK(NCTRK(J)+2,J)  should  always  be  0.0  and  1.0,  respectively.  The 
subscript  I  runs  from  1  to  NCTRK(J)+2.  There  are  no  default 

values. 

EPSK(I,J)  Parameter  specifying  the  degree  of  packing  about  the  point 

XAK(I,J).  It  can  be  thought  of  as  the  length  of  the  interval  about 
XAK(I,J)  in  which  points  a-j  packed.  Smaller  values  of  EPSK 
result  in  tighter  packing.  If  EPSK(I.J),  corresponding  to  the 
location  XAK(I,J),  is  greater  than  the  distance  to  the  location 
XAK(I-l.J)  or  XAK(I+1,J),  no  packing  is  done  about  XAK(1,J). 
Therefore,  since  XAK(I,J)  can  at  most  be  1.0,  an  easy  way  to  not 
pack  at  XAK(I,J)  is  to  set  EPSK(I,J)  =  1.0.  There  are  no  default 
values. 


NAMELIST  VORTEX 

The  parameters  specified  in  this  namelist  define  the  forced  vortex  at  the  initial 
value  surface. 

RCORE  Radius  of  the  core  of  the  forced  vortex. 

QTEDGE  Strength  of  the  forced  vortex  specified  as  the  velocity  at  the  edge 

of  the  core.  Within  the  core  the  vortex  strength  is 

Q=2*QTEDGE/RCORE. 

ZZERO  Location  of  the  center  of  the  forced  vortex  in  the  transverse 

initial  value  surface.  ZZERO  is  a  complex  number,  so  the  y,z 
location  of  the  vortex  is  specified  in  the  form  (y,z). 
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PEPSIG  (CORNER  VORTEX  VERSION) 

ICOEF  OPTIONS 

ICOEF  is  a  6x20  array  read  in  namelist  FLUIDS  that  contains  various  switches, 
limits,  etc.  Some  of  these  apply  to  both  the  potential  flow  (MODE=2)  and  viscous 
flow  (MODE=l)  calculations,  some  to  only  one  or  the  other,  and  some  to  features  in 
the  code  that  are  still  under  development.  For  the  average  user,  many,  if  not  most, 
of  these  options  will  have  little  significance.  Those  most  likely  to  be  needed  by  the 
typical  user  are  marked  with  an  asterisk.  Except  where  noted,  setting  the  value  to 
0  turns  the  option  off.  Except  where  noted,  the  default  values  are  all  zero. 


ICOEF(l.n) 


n  Description 

1  1  Print  initial  station  values  before  starting  procedure. 

2  1  Momentum  integral  boundary  layer  calculation  for  straight 

pipe  flow,  with  the  correct  area  variation,  as  a  basis  for 
computing  the  turbulence  mixing  length. 

The  default  value  is  1. 

3  Maximum  number  of  iterations  in  ADI  for  scalar  potential  solution 

and  for  Poisson  pressure  equation  solution.  The  default  value  is  200. 

4  Maximum  number  of  iterations  in  ADI3D  for  potential  flow  solution. 

The  default  value  is  200. 

5  Maximum  number  of  iterations  in  PRIMRY  for  computation  of  one¬ 

dimensional  viscous  pressure  gradient  correction.  The  default 
value  is  10. 

6  0  Linear  interpolation  on  the  potential  flow  file. 

1  Potential  flow  pressure  coefficient  set  equal  to  0. 

2  Quadratic  interpolation  on  the  potential  flow  file. 

The  default  value  is  1. 

7  Not  used. 

8  1  Print  the  Jacobian  matrix  at  each  g» id  point. 

9  0  Specify  the  centerline  shape  (array  PCL)  by  the  distribution  of 

curvature  and  arc  length. 

1  Specify  the  centerline  shape  (array  PCL)  by  Cartesian 
coordinate  location. 

2  Specify  the  centerline  shape  (array  PCL)  by  Frenet  frame 
parameters. 
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10  1  Print  the  secondary  velocities  in  the  orthogonal  reference 

coordinates  as  computed  during  the  scalar  potential  and 
coupled  vorticity-stream  function  solutions. 

1  1  1  Print  the  final  computed  velocities  in  the  orthogonal 

reference  coordinates. 

12  «£)  Print  iteration  data  in  ADI2X2  during  the  coupled  vorticity- 

stream  function  solution  every  |ICOEF(l,12)|  iterations. 

0  Print  namelist  SCL2X2  in  subroutine  SCALE. 

1  Same  as  0,  plus  print  of  coefficients  of  input  equations  (first 
iteration  only),  plus  iteration  data  every  iteration. 

2  Same  as  1,  plus  print  of  coefficients  of  equations  sent  to  matrix 
inverter  each  sweep,  and  resulting  solution. 

3  Same  as  2,  plus  print  from  matrix  inverter. 

The  default  value  is  -5. 

13  <0  Print  iteration  data  in  ADI  during  solution  of  scalar  potential 

and  Poisson  pressure  equations  every  |ICOEF(l,13)|  iterations. 

0  Print  iteration  data  every  iteration. 

1  Same  as  0,  plus  print  of  namelists  TRMDMP,  DTADI2,  and  DUMP 
in  subroutine  ADI. 

2  Same  as  1,  plus  print  of  coefficients  of  input  equations  (first 
iteration  only). 

3  Same  as  2,  plus  print  of  coefficients  of  equations  sent  to  matrix 
inverter  each  sweep,  and  resulting  solution. 

4  Same  as  3,  plus  print  from  matrix  inverter. 

5  Same  as  1. 

The  default  value  is  -5. 

14  <0  Print  iteration  data  in  ADI  during  solution  of  streamwise 

momentum  every  |ICOEF(  1 , 1 3)|  iterations  (if  ICOEF(5,5)  =  0). 

0  Print  iteration  data  every  iteration. 

1  Print  namelists  TRMDMP,  DTADI2,  and  DUMP  in  subroutine  ADI 
during  solution  of  streamwise  momentum  equation. 

2  Same  as  1,  plus  print  of  coefficients  of  input  equations  (first 
iteration  only). 

3  Same  as  2,  plus  print  of  coefficients  of  equations  sent  to  matrix 
inverter  each  sweep,  and  resulting  solution. 


4 

5 

The 

Same  as  3,  plus  print  from  matrix  inverter. 

Same  as  1. 
default  value  is  -5. 

15 

1 

Print  from  subroutine  FRAME. 

16 

1 

Print  of  initial  profiles  from  subroutine  IPROF,  MXPROF, 
MYPROF. 

or 

17 

2 

Print  coefficients  of  Laplace's  equation  during  potential 
solution. 

flow 

18 

1 

Print  f  'lysical  coordinates,  computational  coordinates, 
differs weights. 

and 
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19  1  Print  of  data  read  from  potential  flow  file. 

20  For  KSTART  =  1  or  2  option,  output  is  printed  during  starting 
procedure  every  ICOEF(l,20)'th  iteration  or  sub-step. 


ICOEF(2,n) 

n  Description 

1  Maximum  number  of  iterations  in  ADI2X2  for  coupled  vorticity- 

stream  function  solution.  The  default  value  is  200. 

2,  3  For  ICOEF(2,3)  greater  than  or  equal  to  zero,  computed  profiles  are 

printed  every  ICOEF(2,3)'th  station  starting  at  station  1C0EF(2,2).  If 
ICOEF(2,3)  is  less  than  zero,  profiles  are  printed  at  stations  specified 
by  the  array  IWSTA.  The  defaults  are  0  for  ICOEF(2,2)  and  1  for 
ICOEF(2,3). 

4,  5  For  ICOEF(2,5)  greater  than  or  equal  to  zero,  computed  results  are 

written  into  the  plot  file  every  ICOEF(2,5)'th  station  starting  at 
station  ICOEF(2,4).  If  ICOEF(2,5)  is  less  than  zero,  results  are  written 
into  the  plot  file  at  stations  specified  by  the  array  IPSTA.  The 
defaults  are  0  for  ICOEF(2,4)  and  1  for  ICOEF(2,5). 

6  If  ICOEF(2,6)  =  1,  the  reference  Cartesian  velocity  components  u,  w, 

and  v  are  put  into  positions  4,  5,  and  6  in  the  plot  file  (overwriting 
KPLT(4),  KPLT(5),  and  KPLT(6)).  These  values  are  useful  for 
particle  tracing  routines. 

7-12  Not  used. 

13  1  Print  non-convergence  message  in  ADI. 

14  1  Non-dimensionalize  all  printed  velocities  by  the  average 

streamwise  velocity  at  the  initial  station. 

15  Used  by  the  code  to  keep  track  of  the  number  of  points  in  separated 
flow  regions. 

16  Grid  point  index  in  the  circumferential  direction  specifying  the 
location  of  the  reference  pressure  when  using  the  external  flow 
option  (IDUCT  =  0). 

17  Grid  point  index  in  the  radial  direction  specifying  the  location  of 
the  reference  pressure  when  using  the  external  flow  option  (IDUCT 
=  0). 

18  1  Print  KGRID  output. 

2  Print  KGRID  output,  plus  KGRID  iteration  data. 

19  1  Print  cross-plane  grid  point  coordinates  when  ICOEF(2,20)  =  1. 
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20  1  Print  geometry  parameters  computed  during  a  MODE  =  3 

calculation  (geometry  set-up  for  potential  flow). 

ICOEF(3,n) 


n  Description 

1-2  Not  used. 

3  Results  are  printed  at  every  (ICOEF(3,3)+l)'th  grid  point  in  the  IY 

(circumferential)  direction. 

4  Results  are  printed  at  every  (ICOEF(3,4)+l)'th  grid  point  in  the  JZ 

(radial)  direction. 

5  Results  are  stored  in  the  plot  file  at  every  (ICOEF(3,5)+l)’th  grid 

point  in  the  IY  (circumferential)  direction. 

6  Results  are  stored  in  the  plot  file  at  every  (ICOEF(3,6)+l)’th  grid 

point  in  the  JZ  (radial)  direction. 

7  0  Ignore  phi-velocity  terms  when  computing  transverse 

velocity  gradients  as  part  of  the  computation  of  transverse 
pressure  gradients  in  the  Poisson  pressure  equation,  and  in 
the  computation  of  the  vorticity  vector. 

1  Include  the  phi-velocity  terms. 


The 

default  value  is  1. 

8-9 

Not 

used. 

10 

0  Potential  flow  pressure  coefficient  from  potential  flow  file. 

1  Potential  flow  velocity  from  potential  flow  file. 

To  get  potential  flow  pressure  coefficient  and  velocity  from 
potential  flow  file.  ICOEF(3,10)  is  automatically  set  equal  to 

ICOEF(5,l  1)  if  ICOEF(5,ll)  =  0  or  1.  If  ICOEF(5,ll;  =  2,  ICOEF(3,10)  can 
be  0,  1,  or  2. 

1 1 

Not 

used. 

12 

1 

Smooth  the  potential  flow  pressure  coefficient  (or  velocity  if 
ICOEF(5,ll)  =  1)  in  the  streamwise  direction. 

13 

1 

Print  normalized  arc-lengths  when  ICOEF(2,20)  =  1. 

14 

1 

Print  absolute  Cartesian  coordinates  when  ICOEF(2,20)  =  1. 

15 

1 

Print  elements  of  Jacobian  grid  transformation  matrix  when 
ICOEF(2,20)  =  1. 

16-20 

Not 

used. 
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ICOEF(4,n) 

n  Description 


1  Maximum  pseudo-time  step  size  used  in  ADI  for  solution  of  the  scalar 

potential,  energy,  and  swirl  equations  is  divided  by  the  factor 

2**ICOEF(4,l). 

2  Minimum  pseudo-time  step  size  used  in  ADI  for  solution  of  the  scalar 

potential,  energy,  and  swirl  equations  is  divided  by  the  factor 

2**ICOEF(4,2). 

3  0  BETA  =  1.0  (backward  streamwise  differencing)  in  ADI. 

1  BETA  =  0.5  (Crank-Nicholson  streamwise  differencing)  in  ADI. 


4  0  Use  local  minimum  time  step  in  ADI. 

1  Use  local  maximum  time  step  in  ADI. 

5  Maximum  pseudo-time  step  size  used  in  ADI2X2  for  solution  of  the 
coupled  vorticity-stream  function  equations  is  divided  by  the  factor 
2**ICOEF(4,5). 

6  Minimum  pseudo-time  step  size  used  in  ADI2X2  for  solution  of  the 
coupled  vorticity-stream  function  equations  is  divided  by  the  factor 
2**ICOEF(4,6). 

7  Maximum  pseudo-time  step  size  computed  by  SCALE  for  the  source 
term  in  the  solution  of  the  coupled  vorticity-stream  function 
equations  is  divided  by  the  factor  2**ICOEF(4,7). 

8  0  BETA  =  1.0  (backward  streamwise  differencing)  in  ADI2X2. 

1  BETA  =  0.5  (Crank-Nicholson  streamwise  differencing)  in 

ADI2X2. 

9  Maximum  pseudo-time  step  size  used  in  ADI3D  for  solution  of  the 

potential  flow  equation  is  divided  by  the  factor  2**ICOEF(4,9). 

10  Minimum  pseudo-time  step  size  used  in  ADI3D  for  solution  of  the 

potential  flow  equation  is  divided  by  the  factor  2**ICOEF(4,10). 

11  - 1  Use  local  mean  time  step  in  ADI3D. 

0  Use  local  minimum  time  step  in  ADI3D. 

1  Use  local  maximum  time  step  in  ADI3D. 

The  default  value  is  -1. 

12  0  BETA  =  1.0  (backward  streamwise  differencing)  in  ADI3D. 

1  BETA  =  0.5  (Crank-Nicholson  streamwise  differencing)  in 

ADI3D. 


13 


Convergence  criteria  in  ADI  for  scalar  potential,  Poisson  pressure, 
energy,  and  swirl  equations  is  multiplied  by  the  factor 
10**ICOEF(4,13). 
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14 

15 

16 

17 


18 


19 


Convergence  criteria  in  ADI2X2  for  the  first  equation  (vorticity)  is 
multiplied  by  the  factor  10**ICOEF(4,14). 

Convergence  criteria  in  ADI2X2  for  the  second  equation  (stream 
function)  is  multiplied  by  the  factor  10**ICOEF(4,15). 

Convergence  criteria  in  ADI3D  is  multiplied  by  the  factor 
10**ICOEF(4,16). 

Number  of  time  step  cycles  used  in  ADI3D  for  iterative  solution  of 
the  potential  flow  equation.  The  maximum  number  is  3.  The  default 
value  is  3. 

1  Base  the  pseudo-time  step  for  the  first  cycle  in  ADI3D  on  a 
combination  of  the  magnitudes  of  the  finite-difference 
operators  in  the  circumferential,  radial,  and  streamwise 
directions. 

2  Base  the  time  step  on  a  combination  of  the  circumferential  and 

radial  directions. 

3  Base  the  time  step  on  a  combination  of  the  radial  and 
streamwise  directions. 

4  Base  the  time  step  on  a  combination  of  the  circumferential  and 

streamwise  directions. 

5  Base  the  time  step  on  the  circumferential  direction. 

6  Base  the  time  step  on  the  radial  direction. 

7  Base  the  time  step  on  the  streamwise  direction. 

The  default  value  is  5. 

Same  as  ICOEF(4,l8),  except  for  second  cycle.  The  default  value  is  6. 


20  Same  as  ICOEF(4,18),  except  for  third  cycle.  The  default  value  is  7. 


ICOEF(5,n) 


n 


Description 


1 

2 


1  Set  the  rotational  part  of  the  cross-flow  velocity  (computed 

during  the  coupled  vorticity-stream  function  solution)  to  zero. 

- 1  Solve  the  Poisson  pressure  equation,  but  get  the  density  from 

the  potential  flow  pressure. 

0  Solve  the  Poisson  pressure  equation,  and  get  the  density  from 

the  viscous  flow  pressure. 

1  Skip  solving  the  Poisson  pressure  equation,  and  get  the  density 

from  the  potential  flow  pressure. 

The  default  value  is  -1. 


3 


Used  by  the  program  in  subroutine  COEFVS.  Do  not  change. 
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4 

5 


6 


0  First  order  wall  vorticity  boundary  condition. 

1  Second  order  wall  vorticity  boundary  condition. 

0  Iterate  the  primary  momentum  equation  to  convergence  each 
marching  step  (mainly  important  in  separated  flow  regions). 
If  used,  convergence  criteria  (controlled  by  ICOEF(5,8))  should 
be  tightened. 

1  No  iteration.  The  default  value  is  1. 

Maximum  pseudo-time  step  size  used  in  ADI  for  solution  of  the 
primary  momentum  equation  is  divided  by  the  factor  2**ICOEF(5,6). 
This  only  applies  if  iteration  is  used  (ICOEF(5,5)  =  0). 


7  Minimum  pseudo-time  step  size  used  in  ADI  for  solution  of  the 

primary  momentum  equation  is  divided  by  the  factor  2**ICOEF(5,7). 
This  only  applies  if  iteration  is  used  (ICOEF(5,5)  =  0). 


8 


9 


10 


11 


12 

13 


14 


Convergence  criteria  in  ADI  for  primary  momentum  equation  is 
multiplied  by  the  factor  10**ICOEF(5,8).  This  only  applies  if 
iteration  is  used  (ICOEF(5,5)  =  0). 

<0  if  estimated  initial  velocity  profiles  are  given  as  velocity 
components  in  the  body-fitted  computational  coordinate  system 
(instead  of  the  centerline  coordinate  system). 

0  No  slip  of  secondary  velocities. 

1  Slip  of  secondary  velocities. 

The  default  value  is  0  for  IBETA  «  1,  and  1  for  IBETA  =  0. 

0  Potential  flow  solution  from  a  pressure  file  (or  to  store  a 
pressure  file  if  MODE  =  2). 

1  Potential  flow  solution  from  a  velocity  file  (or  to  store  a 
velocity  file  if  MODE  =  2). 

2  Potential  flow  solution  from  a  file  containing  both  pressure 
coefficient  and  velocity  (but  this  is  not  currently  available 
through  PFLOW). 

1  Include  rotationality  of  inviscid  velocity  in  COEFVS. 

0  Ignore  phi-velocity  terms  when  computing  streamwise 

velocity  gradients  as  part  of  the  computation  of  transverse 
pressure  gradients  in  the  Poisson  pressure  equation. 

1  Include  the  phi-velocity  terms. 

Maximum  pseudo-time  step  size  used  in  ADI  for  solution  of  the 
Poisson  pressure  equation  is  divided  by  the  factor  2**ICOEF(5,14). 


15  Minimum  pseudo-time  step  size  used  in  ADI  for  solution  of  the 

Poisson  pressure  equation  is  divided  by  the  factor  2**ICOEF(5,15). 
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16 


17 


18 

19 

20 


<0  Print  iteration  data  in  ADI  during  solution  of  energy  equation 
every  |ICOEF(5,16)|  iterations. 

0  Print  iteration  data  every  iteration. 

1  Same  as  0,  plus  print  of  namelists  TRMDMP,  DTADI2,  and  DUMP 
in  subroutine  ADI. 

2  Same  as  1,  plus  print  of  coefficients  of  input  equations  (first 
iteration  only). 

3  Same  as  2,  plus  print  of  coefficients  of  equations  sent  to  matrix 
inverter  each  sweep,  and  resulting  solution. 

4  Same  as  3,  plus  print  from  matrix  inverter. 

5  Same  as  1. 

<f)  Print  iteration  data  in  ADI  during  solution  of  swirl  equation 
every  |ICOEF(5,17)|  iterations. 

0  Print  iteration  data  every  iteration. 

1  Same  as  0,  plus  print  of  namelists  TRMDMP,  DTADI2,  and  DUMP 
in  subroutine  ADI. 

2  Same  as  1,  plus  print  of  coefficients  of  input  equations  (first 
iteration  only). 

3  Same  as  2,  plus  print  of  coefficients  of  equations  sent  to  matrix 
inverter  each  sweep,  and  resulting  solution. 

4  Same  as  3,  plus  print  from  matrix  inverter. 

5  Same  as  1. 

Not  used. 

1  Update  streamwise  velocity  and  turbulent  viscosity  coefficient 
when  using  KSTART  -  1  option. 

Not  used. 


ICOEF(6,n) 


None  of  these  options  are  currently  being  used. 
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PEPSIG  (CORNER  VORTEX  VERSION) 

"Z"  ARRAY  PRINTOUT  AVAILABLE 
WITH  KPRT  OPTION 

This  table  lists  the  variables  that  may  be  printed  and/or  written  into  the 
plot  file  using  the  KPRT  and  KPLT  arrays.  Variable  numbers  1  through  70  are 
stored  in  the  Z  array.  Those  above  70  are  additional  variables  that  may  be  printed 
and/or  written  into  die  plot  file.  In  this  table,  "n"  refers  to  the  upstream  station, 
and  Mn+1"  refers  to  the  station  most  recently  computed. 


No. 

VARIABLE 

MNEMONIC 

NAME 

DESCRIPTION 

1 

u 

NUN 

Streamwise  velocity  at  n 

2 

u 

NU 

Streamwise  velocity  at  n+1 

3 

V 

NVN 

Velocity  at  n  in  relative  Cartesian  y-direction 

4 

V 

NV 

Velocity  at  n+1  in  relative  Cartesian  y-direction 

5 

w 

NWN 

Velocity  at  n  in  relative  Cartesian  z-direction 

6 

w 

NW 

Velocity  at  n+1  in  relative  Cartesian  z-direction 

7 

p 

NRHON 

Static  density  at  n 

8 

p 

NRHO 

Static  density  at  n+1 

9 

Q 

NVORN 

Streamwise  vorticity  at  n 

10 

a 

NVOR 

Streamwise  vorticity  at  n+1 

11 

cp 

NCPIN 

Inviscid  pressure  coefficient  at  n 

12 

cp 

NCPI 

Inviscid  pressure  coefficient  at  n+1 

13 

P 

MUN 

Laminar  viscosity  coefficient  at  n 

14 

P 

MU 

Laminar  viscosity  coefficient  at  n+1 

15 

Pt 

MUTN 

Turbulent  viscosity  coefficient  at  n 

16 

Pt 

MUT 

Turbulent  viscosity  coefficient  at  n+1 

17 

P 

NPRESN 

Static  pressure  at  n  from  Poisson  pressure 
equation 
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18 

P 

NPRES 

Static  pressure  at  n+1  from  Poisson  pressure 
equation 

19 

V 

NVPHN 

Scalar  potential  velocity  in  relative  Cartesian 
y-direction  at  n 

20 

V 

NVPH 

Scalar  potential  velocity  in  relativ  Cartesian 
y-direction  at  n+1 

21 

w 

NWPHN 

Scalar  potential  velocity  in  relative  Cartesian 
z-direction  at  n 

22 

w 

NWPH 

Scalar  potential  velocity  in  relative  Cartesian 
z-direction  at  n+1 

23 

y 

NXYZA 

Reference  Cartesian  coordinate  in  y-direction 

24 

z 

NXYZA+1 

Reference  Cartesian  coordinate  in  z-direction 

27 

X 

NXYZA+1 

Reference  Cartesian  coordinate  in  x-direction 

26-34 

NEI11+ 

Elements  of  Jacobian  matrix  at  n+1 

35-43 

NEN11+ 

Elements  of  Jacobian  matrix  at  n 

44 

2D:D 

NDD 

Dissipation  function 

45 

h 

NHN 

Orthogonal  metric  scale  factor  at  n 

46 

h 

NH 

Orthogonal  metric  scale  factor  at  n+1 

47 

vi 

NVIN 

Inviscid  velocity  at  n  in  relative  Cartesian 
y-direction 

48 

Vi 

NVI 

Inviscid  velocity  at  n+1  in  relative  Cartesian 
y-direction 

49 

wi 

NWIN 

Inviscid  velocity  at  n  in  relative  Cartesian 
z-direction 

50 

wi 

NWI 

Inviscid  velocity  at  n+1  in  relative  Cartesian 
z-direction 

51 

NPX1 

y-component  of  the  transverse  pressure  gradient 

52 

NPX2 

z-component  of  the  transverse  pressure  gradient 

53 


NRAD1 


x-component  of  the  distance  vector  to  the  center 
of  rotation 
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54 

NRAD2 

y-component  of  the  distance  vector  to  the  center 
of  rotation 

55 

NRAD3 

z-component  of  the  distance  vector  to 
of  rotation 

the  center 

56 

V 

NVPSN 

Vector  potential  velocity  in  relative 
y-direction  at  n 

Cartesian 

57 

V 

NVPS 

Vector  potential  velocity  in  relative 
y-direction  at  n+1 

Cartesian 

58 

w 

NWPSN 

Vector  potential  velocity  in  relative 
z-direction  at  n 

Cartesian 

59 

w 

NWPS 

Vector  potential  velocity  in  relative 
z-direction  at  n+1 

Cartesian 

60 

ui 

NUIN 

Inviscid  streamwise  velocity  at  n 

61 

ui 

NUI 

Inviscid  streamwise  velocity  at  n+1 

62 

i 

« 

NLEN 

Turbulence  mixing  length 

63 

V 

NPSI 

Secondary  flow  stream  function 

64 

4> 

NPHI 

Secondary  flow  scalar  potential 

65 

T 

NTEMN 

Static  temperature  at  n 

66 

T 

NTEM 

Static  temperature  at  n+1 

67 

E° 

NEON 

Total  enthalpy  at  n 

68 

E° 

NEO 

Total  enthalpy  at  n+1 

69 

NDPIDX 

Inviscid  streamwise  pressure  gradient 

70 

71 

Mn 

NLOCMA 

Mach  number 

72 

P 

NDSTPR 

Static  pressure 

73 

Pt 

NDTOPR 

Total  pressure 

74 

P 

NDRHO 

Static  density 

75 

CP 

NDCP 

Static  pressure  coefficient 

76 

Tt 

NDTOTM 

Total  temperature 
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PEPSIG  (CORNER  VORTEX  VERSION) 
PARAMETERS 


In  PEPSIG,  the  sizes  of  the  dimensioned  arrays,  and  hence  the  storage  required  for  the 
program,  are  set  using  PARAMETERS.  These  PARAMETERS  themselves  are  set  in  COMDECK 
CPARAM.  Larger  or  smaller  dimensions  can  be  set  for  the  entire  program  simply  by 
changing  the  appropriate  PARAMETERS  in  COMDECK  CPARAM,  and  then  recompiling  the 
program.  The  basic  PARAMETERS  are  defined  as  follows: 

NDYP  -  Maximum  number  of  grid  points  in  the  circumferential  direction  for 
the  viscous  flow  calculation.  Currently  setequal  to  55. 

NDZP  -  Maximum  number  of  grid  points  in  the  radial  direction  for  the 
viscous  flow  calculation.  Currently  set  equal  to  50. 

NDXPA  -  Maximum  number  of  grid  points  in  the  streamwise  direction  for  the 
potential  flow  calculation.  Currently  set  equal  to  50. 

NDYPA  -  Maximum  number  of  grid  points  in  the  circumferential  direction 
for  the  potential  flow  calculation.  Currently  set  equal  to  20. 

NDZPA  -  Maximum  number  of  grid  points  in  the  radial  direction  for  the 
potential  flow  calculation.  Currently  set  equal  to  20. 

MVARP  -  Total  number  of  variables  stored  in  the  Z  array.  Currently  set 
equal  to  70. 


Several  additional  PARAMETERS  are  defined  in  COMDECK  CPARAM  as  functions  of 
those  listed  above.  The  following  PARAMETERS  are  used  in  various  common  blocks, 
DIMENSION  statements,  and  EQUIVALENCE  statements: 


NDP 

= 

Maximum  of  NDYP,  NDZP,  NDYPA,  and  NDZPA. 

NDYNDZ 

= 

NDYP*NDZP. 

NDP2 

= 

NDP**2 

NDPM2 

= 

NDP  -  2 

NDZP1 

= 

NDZP  +  1 

NDPA 

NDZPA. 

= 

Maximum  of  NDYP,  NDZP,  NDXPA,  NDYPA, 

NPFA 


1 8*NDXPA*NDYPA*NDZPA 
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NPF  =  18*NDXPA*NDYPA*NDZPA  -  15*NDYP*NDZP 

-  35*NDYP*NDZP  -18*NDP**2 

-  MVARP*NDYP*NDZP 


The  PARAMETER  NPF  may  require  some  explanation.  The  total  amount  of  storage 
required  in  the  C  array  (common  block  BLKMM)  for  the  potential  flow  calculation  is 
18*NDXPA*NDYPA*NDZPA.  However,  the  total  amount  required  in  common  block  BLKMM 
for  the  viscous  calculation  is  only  15*NDYP*NDZP  (for  array  C),  plus  35*NDYP*NDZP  + 
18‘NDP**2  (for  array  CQVQ1),  plus  MVARP*NDYP*NDZP  (for  array  Z).  The  array 
CPFLOW(NPF)  is  therefore  added  to  common  block  BLKMM  to  make  it  large  enough  for  the 
potential  flow  calculation.  If  BLKMM  is  already  large  enough,  NPF  is  equal  to  1. 

In  addition  to  the  above  PARAMETERS,  the  following  are  used  in  the  BLOCK  DATA 
routine: 


NDC 

=  1 5*NDYP*NDZP 

NZV 

=  MVARP*NDYP*NDZP 

MVARP1 

=  MVARP +  1 

MVARP2 

=  MVARP  +  2 

MVARP3 

=  MVARP  +  3 

MVARP4 

=  MVARP  +  4 

MVARP5 

=  MVARP  +  5 

MVARP6 

=  MVARP  +  6 

MVARP7 

=  MVARP  +  7 

MVARP8 

=  MVARP  +  8 

MVARP9 

=  MVARP  +  9 

MVAR4 

=  MVARP  -  4 

MVAR69 

=  MVARP  -  69 

NDANG 

=  4*NDYPA 

NDRAD 

=  4*NDZPA 

NDC  PI 

=  4*NDYPA*NDZPA 
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RUNNING  A  CASE  WITH  PEPSIG 
(CORNER  VORTEX  VERSION) 


This  input  runstream  is  stored  as  [ LEVY . VOR] VOR . JOB . 


JOB, JN=VOR, MFL=90000Q , T=200 , US=DEFER . 

ACCOUNT, AC=28X162X8, US=LEVY, UPW=PEPSIG. 

COPYD . 

REWIND, DN=$IN . 

ACCESS, DN=OBJECT, PDN=PGB*VOR, ID=SRA . 

ACCESS, DN=UTIL, PDN=LIB*UTILB, ID=SRA, OWN=DE JONG . 

LDR, DN=OB JECT, LIB=UTIL, NA . 

EXIT. 

/EOF 

1  FLOW  WITH  VORTEX 
&RESTRT 
&END 
&FLUIDS 

NS 1-1 , NS2=- 1 , NS 3=1 , NS4=-2 , 

BLD=. 004, 0 . , .02,0., 

VGEN( 1,20) =0.004, 0.000,  0., 

VGEN (1,  19) =0 .010,  0.085,  0  .  , 

REY=1 . E6 , KTURB=1 , IBULEV=0, 

IPLOT=l , AP  =  1 . 07  , 

NEY=4  9,NEZ=49, 

T=0. 038, DTE=. 00811, 

NS=33 , 

ICOEF (1,1) =0, ICOEF ( 1 , 6 ) =1 , ICOEF ( 2 , 1 ) =500 , 

ICOEF (1, 2) =0, ICOEF  (1, 5) =10, 

ICOEF  (1,  12)  =-20,  ICOEF  (1,  13)  =-20,  ICOEF  (1,  14)  =-20, 
ICOEF  (2,3) =-1000, ICOEF  (2, 5)=-1000, 

ICOEF (2, 16) =40, ICOEF (2, 17) =40, 

ICOEF (5, 2) =-l, 

IWSTA=5 , 12,33, 

IPSTA=5 , 12,33, 

ISKPL=1, JSKPL=1,NY1PL=1,NZ1PL=1,  NY2PL=4 5 , NZ2PL=4 4 , 
ISKPR= 1 , JSKPR=1 , NY 1PR= 1 , NZ 1PR=1 ,  NY2PR=4 5 , NZ2PR=4 4 , 
&  END 
&GEOM 
IDUCT=0 , 

VIS= . 004 ,  .9434,  .002,  1., 

NGEOM=2 1 , 

PGEOD  =  1000*0  .  , 

PGEOD (1,1)= . 2, PGEOD (1, 2) =. 116, 

&END 
& VORTEX 

ZZERO=( .023,  .0077)  , 

QTEDGE=0 . 35, 

RCORE=0 . 0023, 

SEND 

STOP 

/EOF 
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WIDTH. HEIGHT, DEL  Y,  DEL  Z:  26  0.20000E+00  0.11600E+00  0.00000E+00  0.0G000E+00 

ADI  **  ITER  I’MAX  DELE  RIISMAX  PTEST  RHS  TEST  LOG  RES 

20  0 . 4  4  672E-03  0.53993E-09  0.98486E-06  0.12087E-01  0.10965E+01  0.37 

30  0 . 44564E-03  0.42333E-09  0.98487E-06  0.94993E-02  0.85966E+00  0.25 
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